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Moduli fields, a natural prediction of any supergravity and superstring-inspired supersymmetry 
theory, may lead to a prolonged period of matter domination in the early Universe. This can be 
observationally viable provided the moduli decay early enough to avoid harming nucleosynthesis. If 
primordial black holes form, they would be expected to do so before or during this matter dominated 
era. We examine the extent to which the standard primordial black hole constraints are weakened 
in such a cosmology. Permitted mass fractions of black holes at formation are of order 10^*, rather 
than the usual 10"^" or so. If the black holes form from density perturbations with a power-law 
spectrum, its spectral index is limited to n < 1.3, rather than the n < 1.25 obtained in the standard 
cosmology. 
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I. INTRODUCTION 

Although a substantial amount of work has been car- 
ried out on the assumption of a 'standard' cosmology, 
in which the Universe proceeds from an early period of 
inflation through reheating to radiation domination and 
finally to matter domination in the recent past, there is 
no direct evidence supporting this picture until the rela- 
tively late epoch at which nucleosynthesis occurs. Re- 
cently, this standard picture has been questioned and 
some alternative cosmologies discussed. An example is 
'thermal inflation' ||l|,|^, a short second period of inflation 
at lower energy scales which does not generate interesting 
density perturbations, but which may resolve additional 
relic density problems not solved by the original infla- 
tionary period. 

Several cosmological constraints are sensitive to what- 
ever assumption is made for the entire cosmological evo- 
lution; axion cosmology is one such situation and 
the constraints originating from primordial black holes 
(PBHs) is another. Recently, the latter was reinves- 
tigated for cosmologies with thermal inflation, showing 
that the standard constraints on the formation den- 
sity of PBHs would weaken quite markedly Q . 

Another possible modification to the standard cosmol- 
ogy is the addition of a prolonged period of matter dom- 
ination, induced by a slow-decaying massive particle. In 

= I supergravity models (b), supersymmetry (SUSY) 
is broken in some hidden sector and the gravitational 
strength force plays the role of messenger by transmit- 
ting SUSY breaking down to the visible sector. In these 
models there often exist scalar fields with masses of the 
order of the weak scale and gravitational strength cou- 
pling to the ordinary matter. If at early epochs one of 
these fields is sitting far from the minimum of its po- 
tential with an amplitude of order of the Planck scale. 



the coherent oscillations about the minimum will eventu- 
ally dominate the energy density of the universe. These 
fields will then behave like nonrelativistic matter, and 
decay at very late times. The presence of these slow- 
decaying massive particles is predicted not only in some 
specific classes of supergravity models, but in almost all 
theories in which supersymmetry is broken at an inter- 
mediate scale. In string models, massless fields exist in 
all known string ground states and parametrize the con- 
tinuous ground state degeneracies characteristic of su- 
persymmetric theories. These fields are massless to all 
orders in perturbation theory, and get their mass, of or- 
der a TeV, from the same non-perturbative mechanism 
which breaks SUSY. Being coupled to the ordinary mat- 
ter only by gravitational strength couplings, a long life- 
time results. Possible examples are the dilaton of string 
theory and the massless gauge singlets of string com- 
pactifications, and they go generically under the name 
of moduli. Under natural assumptions on the couplings, 
one finds that the reheating temperature after the mod- 
uli decay is too low to allow standard nucleosynthesis. 
Therefore, moduli are generally far too good at giving a 
period of matter domination, lasting beyond the epoch 
of nucleosynthesis and destroying this crucial success of 
the standard cosmology (7|j|]. 

Many attempts have been made to resolve this cosmo- 
logical moduli problem However, all of them require 
new phenomena to occur on the cosmological side as well 
as in the theory of supersymmetry breaking. It is not 
our intention, in this paper, to propose another solution 
to the moduli problem, but rather to study the extent 
to which the modular cosmology may affect the stan- 
dard primordial black hole constraints. We will there- 
fore assume that the moduli are somewhat more massive 
than the supersymmetry scale set by the gravitino mass 
"^3/2 = (10^ — 10'^) GeV, and that their decay can be 
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just early enough. This assumption seems reasonable in 
view of the recent developments in the context of the 
field-theory limit of superstrings where it has been 
shown that moduli masses as large as 10'' 7713/2 can be 
achieved without incurring excessive unnaturalness. The 
cosmological sequence in such a model is as follows. At a 
high temperature the moduli come to dominate and the 
Universe begins an epoch of matter-domination. Dur- 
ing this period, the radiation field actually cools to some 
way below the nucleosynthesis scale (about 10~^GeV), 
but the moduli decay while the energy density is still high 
enough to permit thermalization slightly above the nu- 
cleosynthesis temperature. In this picture, baryogenesis 
must be caused by the decay of the moduli, rather than 
at the electro- weak transition [llll. 



II. THE MODULI-DOMINATED EPOCH 

In hidden-sector models, supersymmetry breaking is 
conveyed to the low-energy visible sector through Planck 
scale suppressed interactions. In non-rcnormalizable 
hidden-sector models, supersymmetry vanishes in the 
limit mpi —^ 00, mpi being the Planck mass. Since the 
potential for a generic moduli field (j) is generated through 
the same physics associated to supersymmetry breaking, 
its potential takes the form 



V{^)=ml/^A4iVM/Mpi) 



(1) 



where Mpi — mpi/y/Sn is the reduced Planck mass and 
™3/2 ~ 1 TeV is the gravitino mass. The potential for 
this dangerous direction vanishes in the flat-space limit 
since TO3/2 ^ in that limit. As mentioned in the intro- 
duction, excitations around the zero-temperature mini- 
mum (pQ of the potential have a mass — 0(10^) ^3/2- 

Moduli fields are expected to be initially shifted from 
their zcro-tcmperature minimum due to the effect of ther- 
mal fluctuations or of quantum fluctuations during in- 
flation Another source of the shift might be the 
fact that the moduli couplings to the inflaton generally 
modify, during inflation, the properties of the effective 
potential. Moduli usually acquire a mass squared of the 
order of H^, where H ~ 10^'^ GeV is the Hubble param- 
eter during the inflationary stage, and the value of the 
minimum of the potential may be shifted [ p^ . The shift 
produced by such effects may be as large as mp\. 

Although the form of the potential is not known, for 
our purposes one may just consider oscillations around 
the minimum with initial amplitude and take U(0) ~ 
m'^cj)'^ /2. When the Hubble parameter H reaches a value 
H ^ m^, the scalar fleld starts oscillating coherently 
around the minimum of the potential. This happens 
when the temperature of the universe is Ti ~ ^m^mp\ 
(in the case in which the universe is radiation dominated 
at that epoch). 

The initial energy stored in the oscillations pi ~ 
redshifts like matter and can eventually dominate the en- 



ergy density. When it does so depends on <j)i. If (^i is of 
order mpi, then the moduli dominate immediately, while 
if (j>i is smaller, radiation domination will continue for a 
while before the moduli come to dominate, or, in extreme 
cases, the moduli may decay before they dominate the en- 
ergy density. In hidden-sector models, moduli couple to 
other fields only through Planck suppressed interactions. 
Examples of such fields are the dilaton and the compact- 
ification moduli of string theory or, in general, for any 
gauge singlet field responsible for SUSY breaking. There 
are several types of Planck suppressed couplings the mod- 
uli might have with ordinary matter, but all of them lead 
to the same estimate of the decay width 
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The condition for the moduli to dominate the energy den- 
sity of Universe when they decay (which will be at the 
epoch H ~ r^) is that their initial value satisfies 

At the decay time the radiation fiuid has temperature 
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-1/6 -2/3 



(4) 



The decay products of the moduli will thcrmalize, reheat- 
ing the universe up to a temperature 



3/2 

tpi 



Notice that the reheating temperature is independent of 
0i, provided that the universe is dominated by the moduli 
energy density when decays start. 

The decay products of (f) will destroy the ''He and D 
nuclei, and thus successful nucleosynthesis predictions, 
unless Tioh is larger than about 1 MeV. If the mod- 
uli field has mass 10^ GeV, Treh is well below the en- 
ergy scale of nucleosynthesis, but if instead one assumes 
~ 10'' GeV, then the reheat temperature becomes 
comparable and it may be possible to thermalize to a 
high enough temperature for standard nucleosynthesis to 
proceed. In the case ~ mp\ where the moduli domi- 
nate as soon as they being to oscillate, this corresponds 
to an expansion of the Universe during matter domina- 
tion by a factor of around (mpi/m^)''/'' '-^ 10^°, a very 
prolonged period indeed. 



III. PRIMORDIAL BLACK HOLE 
CONSTRAINTS 

A. Formation density constraints 

In a radiation-dominated Universe at temperature T, 
the horizon mass is given roughly by 
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PBHs of a given mass are expected to form around the 
time when that mass equals the horizon mass; production 
of smaller black holes is suppressed as pressure prevents 
the collapse of any density perturbation. In a matter- 
dominated Universe, formation may occur on scales be- 
low the horizon mass, as we discuss later. 

The lifetime of the black hole can be parametrized as 
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9 X 10 



-27 



f{M) 



M 



sec , 
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where f{M) depends on the number of particle species 
which can be emitted and is normalized to 1 for holes 
which emit only massless particles. A black hole of ini- 
tial mass around 5 x 10^'' g would be evaporating at the 
present epoch, while masses around 10^'' g would be evap- 
orating at nucleosynthesis. Those lighter black holes may 
form early on in the period of moduli domination, or even 
before it if its onset is delayed. 

We denote the fraction of the density of the Universe 
in black holes of a given mass as /3, with (3i denoting the 
initial density at formation. The ratio of the PBH density 
to the density in other forms is denoted a = /3/(l — 

The various limits which can be placed on the PBH 
density are well known ; we shall use the compilation 
given in Ref. There are a range of constraints from 
effects of evaporation, while for more massive black holes, 
M 10^^ g, the only limit comes from their contribution 
to the present density parameter. An additional, less 
secure constraint arises if one assumes that evaporation 
leaves behind a Planck mass relic |15|. 

All these constraints are expressed as limits on the frac- 
tion of the mass of the Universe in black holes at the 
present or at the time of evaporation. To constrain the 
initial mass fraction, one needs to assume a form for the 
entire cosmology back to the formation epoch, given by 
Eq. (|^). Fig. [l] shows the result of carrying this out for 
the standard cosmology, where the Universe was radia- 
tion dominated until very recently ||]. We see that the 
constraints are extremely tight; typically only something 
like 10"^" of the mass of the Universe is permitted to 
form black holes in the standard cosmology. 

We now turn to our main purpose, examining the 
change in the constraints induced by a period of moduli 
domination. In order to attain a reheating temperature 
of Tioh ~ 10"'^ GeV, so that nucleosynthesis can proceed, 
we require ~ 2 x 10'' GeV. In the extreme case of 

^ mpi, the moduli begin to oscillate at temperature 
Tmd 2 X 10" GeV, since 



Pi = m^mpi = —g^ 



MD ■ 



(8) 



Here ~ 250 is the number of degrees of freedom 

in the minimal supersymmetric standard model. In the 
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FIG. 1. The tightest limits on Oi, for the standard cosmol- 
ogy. The mass is in grams. The relic constraint is shown as a 
dotted Une, emphasizing that it is not compulsory. The solid 
lines, from left to right, represent nil production at nucleosyn- 
thesis, deuterium destruction, He-4 spallation, entropy pro- 
duction, gamma ray production, and total density in PBHs; 
see Ref. [pll for details. 



following, we will consider two scenarios: one in which 
^ mpi and moduli domination begins immediately, 
and an intermediate scenario where (pi ~ 10~^ mpi and 
there is a delay before moduli domination commences. 



1. Immediate moduli domination 

From Eq. (||), PBHs in the mass range 2 x 10^ g < 
M < 2 X 10"^^ g are formed during the moduli-dominated 
era. During moduli-domination, the PBHs constitute a 
constant fraction of the total energy density. The time 
of the decay of the moduli, tdoc, and of the reheating of 
the subsequent thermalized fluid, Uch, can be taken as 
the same, giving a PBH density at reheating of 



/ PPBH \ 
V Prad / 



I Ppbh^ 
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/ /^pbh\ 
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(9) 



Here piad and Pmod are the energy densities in radi- 
ation and moduli respectively. Therefore, for PBHs 
formed during moduli domination and surviving beyond 
it (which is all those of interest), we have 



/ PPBH \ _ ^ _ A ^reh 

V Prad / ovap 1 — A Tevap 



(10) 



Considering the duration of the various phases of the 
evolution of the universe 



t, 



cvap 
tpi 



tr 



^dcc 



tj-ch ^Pl ' 



(11) 



using the relation between the formation time and mass 
of a PBH for PBHs formed during radiation domina- 
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tion 



M ^ Mh = —Mpi . 



(12) 



and the variation of the density during moduh domina- 
tion, p cx t~^, leads to 



to. 



T.„h V f Pi M 



tpi V^evap/ VPdcc/ "Ipi 

In order to ehminate pi, Eq. (B) can be rewritten as 
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Mh = 0.2 ■ 
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Finally the resulting expression for iovap can be equated 
with Eq. (^ to give 
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= 8 X 10^21 



\mpi 



3/2 



SO that 
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The gravitational constraints, which require that the 
present-day densities of PBHs and relics not to overdose 
the universe, are 



n 



PBH,oq 



_ / PPBh\ _ A ^reh ^ ^ 
V Prad / cq I ~ Pi T^q 

_ f PPBH \ _ ^Pl A Trch ^ ^ 



(17) 



V Prad / cq Ml- A Toq 

where 'eq' indicates the epoch of matter-radiation equal- 
ity in the Universe's recent past, after which the density 
of PBHs or relics, relative to the critical density, remains 
constant. In the case of PBHs with M > 2 x lO^^^g, 
formed after moduli domination, the requirement that 
the present-day density of PBHs does not overdose the 
universe is obviously the same as in the standard evo- 
lution of the universe. The PBHs formed before moduli 
domination are sufficiently light M < 10^ g that only the 
relic constraint f^rcLcq < 1 applies to them, where: 

n,^,^^=(PEBl) =!!^J^Z^JL<i, (19) 

V Prad / oq M 1 - A T^q TmD 

and in fact it turns out that large initial mass fractions 
of PBHs, A ~ 1) E^i'e allowed. 

The various limits on the initial mass fraction of PBHs 
are illustrated in Fig. |^. 



*This equation is not precisely valid for PBHs which are 
formed during matter domination since their formation may 
be delayed; however, it can be used in this context since the 
delay is negligible compared with the PBH lifetime. We will 
discuss PBH formation during matter domination in more 
detail later. 
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FIG. 2. The tightest limits on the initial mass fraction of 
PBHs, Ui, if moduli domination commences immediately. The 
mass is in grams. The rightmost line, indicating the density 
constraint, continues horizontally until M ~ 10"^* g; PBHs 
more massive than this form after moduli dom ination and 
the standard constraint Qi < 10~^^ -y/ Af/lO^^ g then applies. 
The constraints are the same as in Fig. 1. 



2. Delayed moduli domination 

An initial value 4>i mpi is the most natural, but it 
is not impossible for it to be smaller and this leads to a 
shorter period of moduli domination. As an example, we 
take 0i ~ lO~*m0 so that moduli domination commences 
when the energy stored in the oscillations of the moduli 
field becomes greater than that of the radiation. 



2 X 10"GeV 



Prad 



'PI 



> 1 



(20) 



at temperature Tmd = 2 x 10"^ GeV. 

From Eq. (|), PBHs with M < 2 x 10^^ g are formed in 
the radiation-dominated period before the moduli domi- 
nation commences. Their energy density, relative to that 
in other forms, varies as initially then remains con- 
stant during moduli domination. It then increases as 
during the subsequent radiation domination so that 



/ PPBH \ 

V Piad / cvap 



^cvap 



A Tl T, 



rch 



1 — A Tmd levap 



leading to 



^^2xl0=^^aevap. 

1 - A M 



Similarly for the gravitational constraints 



PBH,cq 



Ppbh 



A Ti Ti-eh 
1 — A ^MD ^e, 



< 1 



Prad / cq ^ ~ Hi -•^ MU J-eq 



(21) 



(22) 



(23) 



VPrad / eq M 1 - A JMD Icq 



4 




log M 

FIG. 3. The tightest hmits on the initial mass fraction of 
PBHs, Qj, if 01 ~ 10~^mpi and moduli domination is delayed. 
The mass is in grams. For Af > 2 x lO^^g the limits are the 
same as when moduli domination commences immediately. 
The constraints are the same as in Fig. 1. 



For PBHs with Af > 2 x 10^^ g, formed after moduli 
domination commences, the gravitational constraint is 
the same as when moduli domination starts immediately 
at T = 2 X 10" GeV. 

The various limits on the initial mass fraction of PBHs 
in this case are illustrated in Fig. ||. 



B. Density perturbation constraints 

PBH formation can be used to constrain the spectral 
index of the density perturbation spectrum |16 l^,^. To 
do this, we assume an initial spectrum which is a power- 
law across the entire range of scales from the PBH scale 
up to the present horizon scale. Constraints can then be 
placed on the spectral index n of those perturbations. To 
do this, we need to consider the variation of aiior{M), the 
mass variance evaluated at horizon crossing 1). During 
matter domination CThor(Af) c>c Af*^^""^/^, whereas during 
radiation domination Chor (M) (X A4^(i-")/4; the different 
scalings arise because during matter-domination the co- 
moving mass density is conserved but during radiation- 
domination it decreases so that mass scales enter the hori- 
zon more quickly 

In the case of PBHs formed during matter domina- 
tion the standard scenario for PBH formation no longer 
holds. It has been shown [Q that perturbation growth 
during coherent scalar field oscillation behaves in exactly 
the same way as in a dust Universe, provided, as here, 
that the oscillations are very rapid compared to other 
timescales in the problem. PBH formation during such 
a matter-dominated period was considered in Ref. [ p^ . 
Because there is no pressure, it is now possible for PBHs 
to form well within the horizon, but in order to do so 
the initial perturbation must be sufficiently spherical as 
gravitational collapse is unstable to aspherical growth. 



The formation rate is given by p9| ] 

/3(A'f) « 2 X IQ-^ a^^/^{M) 



(25) 



For PBHs with A/ > 2 x 10^ g, formed during moduli 
domination, we have 



O-hor(M) = Crhor(Mo) 



M, 



cq 



Mo 



dec 



cq 



(l-n)/6 



(l-«)/4 



(26) 



Af, 



hor 



Mde 



(l-n)/6 



where Mq ~ 10^^ g is the present horizon mass. PBHs are 
formed from rare, relatively large, density fluctuations 
which collapse soon after entering the horizon, so we can 
take Afhor ^ -^^-^ This simplifies to 



-^oAM)^a^oAMo)[j^^j ^^^^^ 



Af 



(l-n)/6 



/Mdo 



(l-n)/12 



(27) 



and since during radiation domination Afn oc T 



O-hor(Af) = 0-hor{Ma) 



/^^(l-«)/6 ^ J. 



oq 



(l-n)/6 



Af 



^dec , 
(l-«)/6 



= ahor(A/o) (^1.4 X lO^'^j , (28) 

for masses Af forming during moduli domination. 

The lightest holes that can form are determined by the 
reheating temperature after the original period of infla- 
tion which is responsible for generating the density per- 
turbations. The minimum mass is then given by Eq. (|^). 
Normally, the tightest constraint on n comes from the 
lightest PBHs. We use the method outlined in ||], but 
using the expressions for cr(Af) and (3{M) given above, 
to obtain the constraints. 

For immediate moduli domination, we find the tight- 
est limit to be n < 1.23 from the deuterium constraint 
evaluated at Af ~ g, although all the constraints 
due to the evaporation of PBHs require n < 1.32. The 
limit from the present-day density of PBHs is tightest at 
Af ^ 5 X 10"'^^ g giving n < 1.30. Relics do not constrain 
n, since even very large initial PBH abundances /3i will 
be diluted away. 

For our example case of delayed moduli domination, 
the most constraining PBHs are formed dviring the 
radiation-dominated era before moduli domination com- 
mences. For them a{M) has a different form 



PBH forming well within the horizon will have a some- 
what smaller mass than this (see e.g. Ref. [Q), but when 
constraining the spectral index the correction is negligible. 
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10 20 30 40 50 



log M 

FIG. 4. The variation of (Jhor(Af) with mass for the three 
scenarios considered; the solid, dotted and dashed hnes rep- 
resenting the standard cosmology, immediate moduli domina- 
tion and delayed moduli domination respectively. The mass 
is in grams, and we take n — 1.3 for illustrative purposes. 

\Mdoc) \MmdJ ^ ' 

which simplifies to 

(Thor(Af) =ahor(Afo) flO'"'^] , (30) 

for our specific parameters. Since the PBHs of interest 
are formed during radiation domination the standard ex- 
pression for /3 applies 

/^(M)..(M)exp(-^^). (31) 

The tightest limit is now n < 1.26 from the deuterium 
constraint evaluated at M ~ 10^° g, with all the con- 
straints due to the evaporation of PBHs require n < 1.28. 
The tightest limit from the present-day density of PBHs 
is n < 1.30 at M ^ 5 x lO^'^g. The relic constraint may 
provide an even tighter limit if the reheat temperature 
after inflation is close to 10^^ GeV. 

In Fig. ^ we illustrate the variation of (Thor (M) in the 
standard cosmology, for immediate moduli domination 
and for our example case of delayed moduli domination. 

The tightest constraint for immediate moduli domina- 
tion is only slightly weaker than in the standard cosmol- 
ogy [||, where the tightest constraint (again deuterium) 
is n < 1.22, with the evaporation constraints all giving 
n < 1.24. There, the relic constraint may give an even 
tighter limit if the reheat temperature after inflation is 
high enough (> lO^'' GeV) to let quite light PBHs form. 
The weakening is only small, since during matter dom- 
ination PBHs form more readily so that to attain any 



particular value of /3i a smaller value of cr(M), and hence 
n, is necessary. This reduces the effect of the larger value 
of Pi allowed due to the period of matter domination, 
and also leads to a larger spread in the limits on n from 
different sources. The tightest constraint is significantly 
weaker for delayed moduli domination, since in this case 
the most constraining PBHs are formed during radiation 
domination so that the main difference from the standard 
scenario is the larger values of f3i allowed. 

IV. CONCLUSIONS 

If there is a prolonged period of matter-domination by 
moduli in the early Universe, it leads to a weakening of 
the constraint on density perturbations from primordial 
black hole formation. It again reminds us of the sensi- 
tivity of this bound to the entire assumed cosmological 
history. If the moduli dominate immediately, the frac- 
tion of the density of the Universe permitted to go into 
PBHs becomes of order 10^*, rather than the 10~^° or so 
which the standard cosmology requires. Delayed moduli 
domination leads to an intermediate constraint on those 
PBHs which form before moduli domination. This weak- 
ening is similar to that found for the case where an 
extra period of inflation at low energies, known as ther- 
mal inflation, is assumed. 

When expressed as a limit on the spectral index of 
a power-law density perturbation spectrum, we obtain 
n< 1.3 for immediate moduli domination, rather than 
n ^ 1.25 as in the standard cosmology. The weaken- 
ing is similar to that from thermal inflation, which also 
led to n < 1.3. Interestingly, the constraint can actu- 
ally be weakest if moduli domination is delayed, because 
PBH formation is harder during radiation domination 
than moduli domination. 

We end by noting that the assumption of gaussian per- 
turbations in the black hole formation calculation has re- 
cently been questioned by Bullock and Primack |Q . As 
shown in Ref. , in the most non-gaussian case found by 
Bullock and Primack the constraint on n can be weak- 
ened further, by up to 0.05. 
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